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I. INTRODUCTION 



A device emitting a coherent atomic beam, similar to a laser radiating coherent photon rays, is called atom laser 
[1-7]. Bose atoms can be prepared in the Bose-condcnsed state by cooling them in a trap [8-10]. The Bose-condensed 
state is believed to be a coherent state [11], which is supported by the interference experiments [12,13]. 

It is worth noting that a Bose-condensed state cannot be a pure coherent state, but it can be only partially coherent. 
This is because a coherent state is not an eigenstate of a Bose system [14]. Another problem is that, if one identifies 
the Bose-condensed state with a state having broken gauge symmetry, then there appear anomalous fluctuations 
of the number of particles, with nonthermodynamic behaviour [15,16]. This anomalous behaviour originates from 
infrared divergences typical of any Bose fluids with broken gauge symmetry [17-19] and is related to the fluctuations 
of phases of field operators [20] . 

However, a system of finite number of particles in a trap does not undergo a genuine phase transition with breaking 
symmetry [21-24], similarly to the case of the Bose gas in restricted geometries [25]. In such cases, the singularities of 
all thermodynamic functions are rounded off and the latter become analytic, varying smoothly. This means that Bose 
condensation of a finite number of particles in a trap is a crossover phenomenon that is not accompanied by breaking 
gauge symmetry. Although breaking the latter can be employed as a technical trick facilitating some calculations, 
for instance, for considering collective excitations of a trapped Bose gas [26-29]. The usage of this trick is limited 
by the cases when it does not lead to anomalies, like those for the fluctuations of the number of particles. The Bose 
condensation in a trap, being a crossover phenomenon, occurs at a crossover temperature that can be defined as a 
temperature at which at least one of thermodynamic characteristics has a maximum. This definition of a crossover 
temperature is typical of crossover phenomena [21-25,30,31]. 

Assume that Bose atoms are trapped and cooled down to experience the Bose condensation, when, at least partly, 
they become coherent. An important problem in realizing an atom laser is how to form a directional beam of atoms. 
To make the output highly directional is, actually, the first condition on a laser [7]. However, at the present time 
there are no mechanisms permitting one to create such well collimated beams of neutral atoms flying out of traps. 
To form an output coupler, one employs short radio-frequency pulses transferring atoms from trapped states to an 
untrapped state. But, when escaping from a trap, atoms fly out more or less in all directions [12,13,32], with an 
inisotropy formed by the gravitational force. This makes the principal difference between the present-day traps and 
optical lasers. "The major difference between the two devices is that the photons from an optical laser generally 
emerge in a well-collimated beam, whereas atoms from the atom laser fly out in all directions" [33] . In order to force 
the escaping atoms to move in one preferable direction, one can employ some mechanical blocking or external laser 
beams (see discussion in Ref. [34]). For example, to provide a directed output beam, an implementation of the laser 
scheme using hollow optical fibers was suggested [35], with a momentum kick provided by two laser beams. 

In this paper, we consider a new mechanism for creating well-collimated beams of neutral atoms. This mechanism 
does not require any mechanical blocking or additional laser operation. Collimation of an atomic beam is achieved only 
by means of magnetic fields of a trap. A particular configuration of the trap magnetic field, leading to a semiconfining 
regime of motion, has been studied earlier [36-38]. Here we generalize the consideration showing that there exists a 
large class of magnetic fields permitting one to create directed beams of atoms for atom lasers. 



When the space variation of magnetic fields is sufficiently smooth, the dynamics of neutral atoms is well described 
in the semiclassical approximation [36,37] for the quantum-mechanical averages of the real-space coordinate, R=< 

R >= {R a }, and of the spin operator, g=< S >= {S a }, where a — x,y,z. Then, for an atom of mass m and 
magnetic moment jio, the equation for the average space variable is 



II. EVOLUTION EQUATIONS 



d 2 R, 



dt 2 



a 




(i) 



and the evolution of the average spin is given by the equation 



d S 
dt 



= TtSxB . 
n 



(2) 



The total magnetic field 
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B=Bi + B 2 

consists of two terms. One is the quadrupole field 

Bi= B[ (r x ~e x +R y e y +XR Z ~e 



(3) 
(4) 



in which 6 a are the unit Cartezian vectors and A is an anisotropy parameter. If one requires that Bi satisfies the 

equation V ■ Bi= 0, then A = —2. But if the field (4) is produced by several different magnetic coils, then one can, 
in general, realize any anisotropy with an arbitrary parameter A. The second term in (3) is a transverse field 



B2= B 2 h (t) , h (t) = h x e x +h y e y , 

in which h a = h a (t) and 

\h\ 2 = h 2 x + h 2 y = l, h = h (t) . 
It is convenient to work with the dimensionless space variable 

r= — = {x,y,z}, R Q = 

Kq tf 1 



Introduce also the characteristic frequencies 



mi? 



UJ 2 



Then equation (1) can be written as 



d 2 r 



A (S x e x +S y e y +XS Z e z ^j . 



This equation is to be complemented by initial conditions 

r* (0) = r*o= W Do, z } , r (0) = r* = {xo, Vo, ^o} 

in which the dot means the time derivative. 

The equation (2) for the average spin can be reduced to the form 



d S 
dt 



uo 2 A S , 



where the matrix A — [A a p] consistes of the elements 

A aj3 = -Ap a , A aa = {a, 13 = 1,2,3) 
A 12 = Xz , A 23 = x + h x , A 31 =y + h y . 

The initial condition for (10) is 

5(0) =s ={S x ,S° y ,S° z }. 



Assume that the following inequalities hold true: 

\ui I <C U) 2 



d h 



dt 



< U) 2 ■ 



Defining an effective frequency u = ui(t) of the transverse field, 

dh 



dt 



we may rewrite the inequalities (13) as 





«1, 


U) 






UJ 2 




LU 2 



< 1 . 



(5) 



(6) 
(7) 

(8) 
(9) 

(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



The existence of these inequalities makes it possible to find approximate solutions to Eqs. (8) and (10) by employing 
the scale separation approach [39-41]. 
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III. SCALE SEPARATION 



The first step in the scale separation approach [39-41] is to classify the functions under consideration into fast and 

slow. Inequalities (15) suggest that the variables r and h are slow, as compared to the fast variables S- It is worth 

noting that equation (8) for r is a differential equation of second order, while that for S is a first-order differential 
equation. Hence, the set of these differential equations does not make what is called the standard form of a dynamical 
system [42], where all equations are to be first-order differential equations. Then one may ask whether it is admissible 

to treat the variable r as slow with respect to S, when not their first derivatives are compared. However, it is not 
difficult to show, by a simple change of notations, as is done in Appendix A, that it is really the case: the variable r 
is low as compared to S- 

With the slow variables treated as quasi-integrals of motion, the solution of Eq. (10) can be found in the following 
way. Solve the eigenproblem 

A bi= a t bi , (16) 

resulting in the eigenvalues 

a h2 = ±ia , a 3 = , a 2 = A\ 2 + A\\ + A 2 zl , 

and in the eigenvectors 

~bi= — ^= [(AlaAas - a t A zl ) ~e x + (A 12 A 31 + a t A 23 ) ~e y + (A 2 12 + af) e . 



fCi 

with the normalization constant Ci given by the equation 

a = (^ 2 -H 2 ) 2 + (^ 2 + K| 2 )(A| 3 + ^ 1 ) . 

It is possible to check that the eigenvectors b » satisfy the orthonormality condition 

b* ■ ~bj= k 3 • (17) 

When the transverse field h is kept constant, then a particular solution of (10) has the form of bi cxp(aiUi 2 t), and 
the general solution is 

-» -» 
S (t) = 2_/fli bi ex P( a i w 20 {h = const) . 

i=l 

Substituting here the dependence of h = h (t) on time and using the notation 

6, (t) =~bi (h (*)) , on(t) = a, (h (t)) , (18) 
we obtain an approximate solution of Eq. (10) in the form 

3 

S (t) = a * bi (t) exp {oti(t)w2t} • (19) 
»=i 

The coefficients a>i can be defined from the initial condition (12), with the use of the orthonormality condition (17), 
which gives 

at =So ■ ~bi (0) . (20) 

A slightly different presentation of an approximate solution to Eq. (10) can be found as follows. Let us look for 

a particular solution of the type bi exp(y>j), keeping in mind that bi and ipi can depend on time. Substituting this 
particular solution into (10) yields 
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bi<pi+ bi = U2Cti bi , (21) 

where the dot means the time derivative. Multiplying (21) by 6 - , we get 

►* _l» 

& = u)2tti— bi - bi , 

from where 

<p(t) = J (w20i- b* -b^j dt . (22) 
Substituting (22) back into (21) gives the equation 

bi = (l>* -~bi^ , (23) 
playing the role of the criterion for the validity of (22). If (23) is valid, then the solution to (10) writes 

3 

S(t) = Y^ <H ~bi (*) exp{v>i(i)} • (24) 

i=l 

Condition (23) imposes a restriction on the behaviour of b % and it may be not valid exactly for an arbitrary field 

h (t). However, here we have the parameter e = |a>/w 2 |, which, accroding to (15), is small, e -C 1. As far as bi 

depends on time through h (t), then bi — > when e — > 0. Therefore, equation (23) becomes asymptotically exact 

when e — > 0. And for small e <C 1, equation (23) is asymptotically valid for an arbitrary field h (t) satisfying (15). 
Thus, the asymptotic form (24) is an approximate solution to equation (10). With the small parameter e < 1, the 
phase (22) is asymptotically equal to 

tp(t) ~ L02a t {t)t (e<l), 

which reduces solution (24) to the form (19). 

In this way expressions (19) and (24) are both approximate solutions to equation (10), provided that the inequality 
e <SC 1 is valid. When e — > 0, these solutions tend to the exact solution corresponding to e — 0. Because of this, such 
solutions are often called asymptotically exact with respect to e — > 0. 

The found solution for the fast variable S has to be substituted into the equation (8) for the slow variable r , 
averaging the right-hand side of this equation over time according to the rule 

</>= Mm - [ f(r,t)dt, (25) 



with the slow variable r kept fixed. This averaging procedure transforms (8) to the equation 

d 2 7 



--F , (26) 



dt 2 

with the averaged force 

F= loI < S x 7 X +S y 7 y +XS Z ~e z > . (27) 
Note that any form of 5", either (19) or (24), can be used, since in both these cases we come to the equality 

<5>=a 3 <T 3 >, (28) 

in which 

_ (x + h° x )S° x + (y + h° y )S° + XzS° _ (x + h x ) 7 X +( y + h y ) l y +Xz e 2 

° 3 l(x + h») 2 + (y + hV) 2 + \ 2 z 2 ]V 2 ' 3 l(x + h x ) 2 + (y + h y ) 2 + \ 2 z 2 }i/ 2 ' 1 ' 

with hP a = h a (0) and a — x,y, z. In this way, for the averaged force (27) we have 

F= uj 2 a 3 < bl ~e x +b\ ~e y +\b z z 7. z > , (30) 
where 63 , with a = x,y, z, are the components of the vector 63= {63}. 
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IV. ROTATING FIELD 



To proceed further, let us concretize the transverse field h (t) in (5). Take this field in the form of the rotating 
field with 

h x = cosuut , h y = sinuit , (31) 
which is employed in some magnetic traps [43,44]. Then, using h x = 1, h® = 0, and the equality 

< xcosojt + y smut >= , 

for the force (30) we find 

- + x)Sl + yS° y + XzS^jx 7, x +y +2A 2 z e z ) 

2[(l + 2x + x 2 +y 2 + A 2 z 2 )(l+a; 2 +y 2 + A 2 z 2 )] 1 /2 ' [ > 

If we take for the initial spin polarization the standard initial condition S® 7^ 0, S® = S® = that is used for 
confining atoms in a trap, then, for 5° < 0, the force (32) really provides confinement, when atoms oscillate in a 
nearly harmonic potential. This type of the confined oscillating motion does not change much if we add to (32) the 
gravitational force. A vertical field gradient supplies the levitating force to support the atom against gravity. The 
combination of the magnetic field and gravity produces a very nearly harmonic confining potential within the trap 
volume in all three dimensions, so that the atom oscillates around an effective equilibrium position [45,46]. The role 
of the gravitational force is discussed in Appendix B. 

We take here a different type of the initial spin polarization corresponding to the initial condition 

S° x = S° y =0, S° z = S^0. (33) 

The choice of initial conditions, as is well known from quantum mechanics, is not prescribed a priori but can be 
realized in any desirable way. Some more details on the possibility of realizing different initial conditions for the spin 
polarization are considered in Appendix C. 

To simplify the evolution equations, it is convenient to measure time in units of u>q , where 

u, 2 = \Su;l\ . (34) 

In what follows, we shall deal with this dimensionless time. To return to the dimensional time, we need to set t — > u>ot. 
In order to take into account the finite size of a trap, we introduce the shape factor 



<p(?) = exp (x 2 + y 2 + e 2 z 2 ) J , (35) 

in which L is the characteristic radius of the device in the radial direction and L/e is the length of the device in the 
axial direction. And let us define the function 

trt) = Pfi f 36 ) 

JK ' [(l + 2a; + a; 2 + y 2 + A 2 z 2 )(l+a; 2 +2/ 2 + A 2 z 2 )] 1 /2 • <• ' 

With these definitions, equation (26), for the initial spin polarization (33), results in the system of equations for 
the components 

^(sgnS)/^, J| = \ (sgnS) fyz , ^ = A 3 (sgnS) fz 2 . (37) 

These equations are invariant with respect to the transformations A — > —A, S — > —S as well as to the transformations 
S — > — S 1 , z — > —z. Therefore, it is sufficient to study only one case for which the sign of XS is fixed, since the change 
\S — ► —XS leads to a picture that is mirror symmetric with respect to the x — y plane. Take, for instance, 

XS>0, \X\=(3. (38) 

Changing AS* by —AS* would invert the whole picture according to the tranformation z — > — z. We may also notice 
that the system of equations (37) possesses the integral of motion 
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d y 

xy— In - = const , (39) 
at x 

which shows that the motion along the x and y axes are similar to each other, so that, under the same initial 
conditions for x and y, the laws of motion x(t) and y(t) coincide. Because of this, we consider in what follows only 
the x component. Thus, taking account of (38), from (37) we have 

£-!'«• $=^ 2 - (40) 

At the initial stage of the process, when | r | -C 1 and f(r) ~ 1, the equations (40) can be solved analytically. 
Then the second equation from (40) writes 

£-flV. (41) 



This can be integrated once yielding 



where 



dz\ 2 2 



dt 3 



= ^ (z 3 - zl + C) , (42) 



3 ,2 



(43) 

Eq. (42) is the Weierstrass equation whose solution 

z(t) = -^V(t-t ) (44) 

is expressed through the Weierstrass function [47]. The time to plays the role of the escape time [36,37] and is defined 
from the initial condition z(0) = zq, which gives 



3_ f 00 dz_ 



J Zn 



The first equation from (40) takes the form 

dx 2 3 

^2- = ^(*-*o)^, (46) 
which is the Lame equation [48] whose degree is defined by the equality n(n + 1) = 3//3 2 , with n > 0, which yields 



-Hv 1 *?- 1 )- < 47 » 

The solution to (46) can be expressed through the Lame functions [48]. 

From the properties of the Weierstrass functions [47] it follows that the axial motion is bounded from below by the 
minimal value 

z min = (z 3 - () 1/3 . (48) 

When t — > to i then 

x{t) ~ |t-t |" 1/2 , a(t) ~ |t-t |" 2 • (49) 

The aspect ratio 
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shows that the atom trajectories are stretched in the axial direction. 

The asymptotic expressions (49) give, of course, only a qualitative understanding of the collimation process, since 
equations (41) and (44) are valid only for the initial stage of atomic motion, when x and z are small. To consider the 
dynamics of atoms for arbitrary times, we have to solve Eqs. (40) numerically. We accomplished such a numerical 
investigation for the parameter (3 = 2 and e = 1. The initial positions of atoms are taken in the center of a trap, with 
varying initial velocities |±o| < 1 and |io| < 1. The atomic trajectories in the x — z plane are presented in Figs. 1 to 
3. The phase portraits for the velocities 

v(t) = % , «,(*) = § (50) 

are shown in Figs. 4 to 6. These figures demonstrate that the bunch of atomic trajectories is essentially squeezed in 
the radial direction forming a well collimated beam. 

V. CONSTANT FIELD 

The limiting case of a transverse field satisfying the second inequality in (15) corresponds to a constant field with 

h x = const , h y = const . (51) 

Then the effective frequency (14) is zero. The matrix A in (10), with the elements (11), does not depend explicitly 
on time. Under r fixed, the solution (19) becomes an exact solution to equation (10). For the average force (30), we 
find 



- 2 (x + h x )S x + (y + hy)S y + \zS° z 
Ul (x + h x )2 + (y + h y )2 + \ 2 zi 



(x + h x ) e x +(y + h y ) e y +\ 2 z e z . (52) 



Note that for the initial spin polarization S x ^ 0, S y = S® = the force (52) does not provide confinement, 
contrary to the average force (32) related to the rotating field (31). However, confinement is not our concern here. 
We wish to find a regime of semiconfinemcnt, when the atoms move predominantly in one direction. 

Let us take for the initial spin polarization the initial condition (33). Then the average force (52) becomes 

F= ^ X(s g nS) Z ( - ± lg |^ + y |*+*| if , (53) 
(a; + h x y + (y + h y y + X 2 z 2 

where the notation (34) is used. For the trap form factor, we accept the same expression (35). But, instead of (36), 
we now have 



(x + h x y + (y + h y y + yzi 



f(r) = - | u - q h , \ 2 .,2 ( 54 ) 



Passing to the dimensionless time measured in units of u> 1 , we obtain as the equation of motion (26) the system of 
equations 



^ = X(sgnS)(x + h x )f Zl ^ = A( hS i..S')(,y -/,,/)/-. 



^=A 3 (sgn5)/z 2 . (55) 

Eqs. (55), though are different form (37), possess the same invariancc property with respect to the transformations 
A — > —A, S — > — S and S — » — S, z — > — z. Thence, we can again accept condition (38). The integral of motion (39) 
exists for (55) only if h x = h y . 

The system of equations (55) also provides the semiconfining regime of motion. For example, at the initial stage, 
when | r | <c 1 and f(r) ~ 1, assuming condition (38), we have 

dX (IZ o n 

d*=? Z - 
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The equation for the y-component is the same as for the variable x, with the change of h x by h y . The equation for z 
coincides with (41), thus, having the same solution (44), demonstrating that the atomic motion is locked from below 
by the minimal value (48). When the time approaches the escape time (45), then instead of (49), we now have 

x(t) ~ln|t-t | , z(i)~ \t-to\~ 2 , 

which shows that, with the transverse constant field, the collimation is even better than with the rotating field. At 
the late stage of the process, when \x/h x \ 3> 1 and \y/h y \ >• 1, the factors (54) and (36) are asymptotically equal, 
and the evolution equations (55) and (37) become equivalent. 

In this way, the semiconfining regime of motion can be realized for different magnetic fields satisfying condition (15), 
under the initial spin polarization (33). This regime is clearly illustrated by numerical calculations for the rotating 
transverse field. The numerical investigation of motion for the constant transverse field will be given in a separate 
paper. The semiconfining regime, we described, can be used for creating well-collimated beams from atom lasers. 
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Appendix A. Standard Form 



Equations (8) and (10) can be written in the standard form, as a system of differential equations of first order, 

d7 -* dp 2 dS , 

with finite functions /i and /2. By means of the notation 

t' = U) 2 t , P = — , £= — , 
U)i UJ 2 

the above differential equations can be presented as 

d7 _ ^ dp'_ dS_ 

-dv~ £P > "^"- e/l ' "^- /2 - 

Condition (15) tells us that £ is a small parameter, |e| <C 1. Then from the latter system of equations it follows 
immediately that r is a slow variable with respect to S- The same concerns the effective velocity P . 



Appendix B. Gravitational Force 

Since an atom has a mass, then, contrary to photon lasers, the existence of gravity may influence the performance 
of atom lasers. The influence of gravity is not essential for the confining regime, when we start with the initial spin 
polarization 5° < 0, 5° = S° = 0. Then, adding to the evolution equation (8) the gravitational force —mg e z , in the 
case of the rotating field, we come to the equations 

d 2x _ l|cO|, ,2- dz_ »2|o0|, ,2„ 9 

dt 2 T xf 1 dt 2 1 xf 1 i? 

describing the motion of atoms inside the trap, where | r | <C 1. We do not write down the equation for y, which is 
similar to that for x. These equations, as is evident, define an oscillatory motion in all three dimensions. The sole 
thing that is changed, when taking account of gravity, is the shift of the equilibrium position on the z-axis from zero 
to z eq = — g(X 2 \S^,\uj 2 Ro)^ 1 . So that the center of harmonic oscillations is shifted from the center of coordinates to 
the point {0, 0,z eq }. Certainly, if the quadrupole field is switched off, the atoms fall down because of gravity. But 
while the quadrupole field is sufficient for confinement, the gravitational force does not change principally the regime 
of motion corresponding to simple harmonic oscillations. Recall that the constant transverse field does not provide 
confinement in any case, so that the gravity is again of no importance. 

In the case of the semiconfining regime, with the initial spin polarization (33), the axis z is an axis of the device, 
but not necessary the vertical axis. Since the device can be oriented arbitrarily, the gravitational force can also be 
directed along different axes. For instance, assume that this force is — mg e x . Then, for the rotating transverse field, 
instead of (40), we have 

d 2 x (3 d 2 z 3 2 

■d¥-r fxz -^ d^ = Pfz > 

where 7 = g/uj 2 R . To estimate the value of 7, let us take the parameters of the magnetic fields as in the quadrupole 
traps with the rotating transverse field [43,44]. Then lu ~ 10 2 - 10 3 s _1 and R ~ 0.1 - 0.5cm. With g w 10 3 cm 
s~ 2 , this gives 7 ~ 10~ 3 — 1. The equation for z is the same as earlier and describes the semiconfining motion along 
the 2-axis. The motion in the radial direction is also scmiconfined, with the downward deviation of an atomic beam 
because of gravity. Thus, the regime of motion does not change principally. 

If the device is oriented so that its z-axis is in the vertical direction, then, instead of (40), we get 

d 2 x /3, d 2 z aZ „ 2 
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Integrating once the equation for z, we have 



(dz_ 
\dt 



2 , 



z - z„ 



^{z-z ) + ( 



The latter equation can be reduced to the Weierstrass form 



eft / 



47> 3 - fl2 P - 53 



with the Weierstrass invariants 



92 



127 



53 



27Zo 



For z, then, we find 



z(t) 



6 



7>(t-t ) • 



When g| < 27g|, then the motion along the z-axis is, as earlier, semiconfined. But if g\ > 27gg, then the motion of 
atoms can become confined depending on initial conditions for z. Recall again, that we are not obliged to align the 
device axis along the gravitational force. Therefore the semiconfining regime of motion can always be realized. 



Appendix C. Initial Conditions 

The evolution equations (1) and (2), we started with, are the equations for the averages of the position and spin 
operators. The average of an operator Q of an observable quantity is the scalar product < Q >= (tp,Qip), where 
tp = ip{t) is the wave function satisfying the time-dependent Schrodinger equation. The evolution equation for such 
an average reads 

d 

ih— < Q >=< [Q,H] > , 

where H is the Hamiltonian (see e.g. [49,50]). The initial state ipo = i/j(0) of the time-dependent Schrodinger equation 
is, in general, an arbitrary function. Consequently, the initial value (iPotQ^o) is a ls° arbitrary. One tells that the 

initial state can be prepared [50,51]. The evolution equation for the average S = S (t) =< S > of the spin operator S 

is equation (2) describing the spin precession from an arbitrary prepared initial value So=S (0), as is discussed, e.g., 
in Ref. [52] . In this way, according to the basic principles of Quantum Mechanics, the initial state can, in general, be 

chosen arbitrarily, resulting in an arbitrary value of 5o- 

It is a different question how this or that particular state can be realized in experiment. One distinguishes two 
limiting cases of preparing initial conditions, adiabatic and nonadiabatic [53,54]. In the first case, the time-dependent 
Schrodinger equation with a Hamiltonian H(t) is ascribed an initial state ipo that is an eigenstate of H(0), i.e., 
H(0)ip = £(0)^0- I n the second case, a chosen initial state tpo is not an eigenstate of H(0). This implies that 
the state ipo could be prepared, for t < 0, as an eigenstate of another Hamiltonian, say Ho, and then, at the time 
t = 0, some external field is suddenly turned on, instantaneously changing the Hamiltonian from Hq to H(0), so that 
Hq ^ H(0). These two types of initial conditions correspond to different physical realizations of an adiabatically slow 
changing field and of a sudden switch on of an external field [53,54]. Both these cases can, in principle, be realized in 
experiment. Of course, when one aims at trapping atoms, the adiabatic motion is necessary [55], with adiabatic initial 
conditions. But this is not mandatory when one's aim is to obtain a specific nonconfined regime. Then one should try 
nonadiabatic initial conditions. Our conditions (33) for the initial spin polarization are such nonadiabatic conditions 
providing the semiconfining regime of motion which can be used for creating well-collimated beams of atom lasers. 
This kind of initial conditions can be realized in practice in different ways. For example, one can confine atoms in a 
trap with a vertical bias field [45], where all atoms become spin-polarized along the z-axis. After this, the bias field 
in the z-direction is suddenly switched off and, at the same time, a transverse bias field, as in experiments [43,44], 
is suddenly switched on. Another possibility could be to prepare spin-polarized atoms in one trap and to suddenly 
load them into another trap with the desired field configuration. We are not going to plunge into technical details of 
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preparing such initial conditions for realizing the scmiconfined motion. This is neither our aim nor our speciality. To 
do this is the challenge for experimentalists. We think that this semiconfined motion can be realized, since there are 
no principal theoretical obstacles for it. 
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Figure Captions 



Fig. 1. The trajectories of atoms at the initial stage of the collimation process, < t < 5, for a trap with L = 100: 
(a) isotropic initial conditions, |x | < 1, |io| < 1; (b) slightly anisotropic initial velocities, |i | < 0.25, \z n \ < 1. 

Fig. 2. Atomic trajectories for the time < t < 20 for two difefrent traps: (a) L = 10; (b) L = 100. 

Fig. 3. Atomic trajectories for long times, < t < 100, for a trap with L = 100. 

Fig. 4. The velocities of atoms in the radial, v(t), and axial, w(t), directions, with initial velocities |±o| < 1 and 
|i | < 1 for two cases: (a) < t < 100, L = 100, (b) < t < 200, L = 10000. 

Fig. 5. Atomic velocities for the initial conditions \xq\ < 0.1 and \zq\ < 0.1 for the following cases: (a) < t < 
100, L = 100; (b) < t < 200, L = 1000. 

Fig. 6. The velocities of atoms for the initial conditions \xq\ < 0.1 and \zq\ < 0.1 for L = 10000, at different stages: 
(a) initial stage, < t < 10; (b) long times, < t < 200. For t > 200, the picture practically does not change. 
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